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o ■ 

We investigate the third-order density perturbation and the one-loop correction to the 
linear power spectrum in the dark-energy cosmological model. Our main interest is to un- 
derstand the dark-energy effect on baryon acoustic oscillations in a quasi-nonlinear regime 
. (k « O.l/i/Mpc). Analytical solutions and simple fitting formulae are presented for the 

dark-energy model with the general time- varying equation of state w(a). It turns out that 
the power spectrum coincides with the approximate result based on the EdS (Einstein de- 
Sitter) model within 1% for k < 0.4/i/Mpc at z = in the WMAP (Wilkinson Microwave 
Anisotropy Probe) 5yr best-fitting cosmological model, which suggests that the cosmological 
i ' dependence is very small. 

O ■ §1. Introduction 

U: 

Revealing the nature of dark energy is fundamentally important not only for 
astrophysics but also for particle physics. Constraints on the dark energy from 
', astronomical observations is very influential for them. Baryon acoustic oscillations 

' (BAO) in the galaxy power spectrum provide a strong constraint on the dark energy 

using its acoustic scale as a standard ruler. Large galaxy surveys such as the Sloan 
Digital Sky Survey and two degree field already provide the constraint and future 
larger surveys are currently planned to detect the BAO more accurately. x )~ 4 ) Hence, 
\Q • an accurate theoretical model of the BAO is crucial, and many authors have been 

investigating the BAO using numerical simulation 5 ) -11 ) and the perturbation theory 
(including the renormalized perturbation theory). 13 )~ 21 ) 

Previously, several authors investigated the third-order density perturbation and 
. ,— i \ derived the one-loop correction to the linear power spectrum in the EdS model. 23 ) -28 ) 

Similarly for the cosmological constant model, Bernardeau (1994) presented the 
third-order perturbation solution (see also Refs. 32)-36)). They found that the 
dependence of the cosmological model on the second- and third-order perturbations 
is very small, if the scale factor in the EdS model is replaced with the linear growth 
factor I*)] However, since the theoretical model of the BAO should archive the sub- 
percent accuracy to provide a strong constraint on the dark energy, it is useful to 
reinvestigate this topic to accurately check the above assumption. In this study, we 
calculate the third-order density perturbation, newly including the dark energy with 
the time-varying equation of state, and derive the one-loop power spectrum analyt- 
ically for the first time. We compare our results with the approximate results based 
on the EdS model in detail, and discuss the effect of the dark energy on the power 



OO 

o 



*' Martel & Freudling (1991) and Scoccimarro et al. (1998) showed that this assumption is valid 
if / = din Di/d In a = fi\{ 2 ■ However, since / » fi^f , the approximation is not accurate. 
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spectrum near the baryon acoustic scale. 

Throughout this paper, we use 6 as the density fluctuation, 6 (= V • v) as the di- 
vergence of the peculiar velocity field, and r = a(t)dt as the conformal time. SIm, 
and fix are the density parameter for the matter, the curvature and the dark energy 
at present. w(a) is the equation of state of the dark energy. The Hubble expansion 
fiua~ z + f2 K a~ 2 + fix exp 3 da' (1 + to (a')) /a' 



rate is H 2 {a) = Hi 



52. Basics 



The equation of motion determines the growth of the density field 5(k,r), and 
velocity field 9(k,r) (= ik ■ v(k,r)) in the Fourier space is 31 ) 



8S(k,r) 



+ 9(k, t)=- J d 3 q a(q, k - q)9(q, r)5(k - q, r), (2-1) 



dr 

T + a(T)H(T)9(k, t) + ^a 2 (T)f2 M (T)H 2 {T)5(k, r) 



9r 



y d 3 q /3(g, fe - g)0(g, r)0(fc - g, r), (2-2) 



with 



(p + g)-p \ (p + g) 2 pg , 00 , 
a(p,g) = 2 — ' P(p>q) = — tti — • ( 2 ' 3 ) 

Equation (|2-ip is the continuity equation, while equation (|2-2p is the Euler equation 
with the Poisson equation. In the linear regime, one can neglect the mode-coupling 
terms on the right-hand sides of Eqs. (|2-ip and (|2-2p . Then the linear solutions are 

Ji(Jfe,a) = Di(a)5i(fc), 0i(fc,a) = -a 2 ff(a) £i(fc). (2-4) 

The linear growth factor Z?i(a) is determined by 

d Dt ( . dhxH\ d Di ( n dlnH 3^ , A D 1 n , n . 
+ 4 + -nr- hnr-— + 3 + -nr-- o^M(a) — = 0, (2-5) 



d 2 In a 2 a \ d In a / d In a a \ d In a 2 / a 

with the initial condition D\(a)/a — > 1 at o — ► 0. In the special case of the flat 
model = 0) with the constant w, the solution is given by the hypergeometric 
function. 37 )' 38 ^ 

The density field is formally expanded up to the third order as 5(k, a) = 5i(k,a)+ 
5 2 (k, a)+53(k, a). We will show the second- and third-order solutions in the following 
sections. 

§3. Second-order solution 

Inserting the linear-order solutions of <!>i and 9\ into the right-hand sides of 
equations (|2-ip and (|2-2p . one can obtain the second-order solution as 



5 2 (k, a) = D 2A (a)A{k) + D 2B (a)B(k), 



(3-1) 
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with 

A(fc) = | | d 3 qa(q,k-q)5 1 (q)5 1 (k-q), (3-2) 

B(fc) = 7/ d3 9 /%> fe - q)Si(q)St(k - q). (3-3) 

The second-order growth factors Dia,b are determined by ordinary differential equa- 
tions with the boundary condition T>2A,B — > a 2 at a — > (see Appendix A). One 
usually approximately use D 2 , instead of D2A,b, in equation (|3-1[) . In order to 
demonstrate the validity of this approximation, we show the relative differences be- 
tween Dia,b and D\ in Fig. [I] for the constant w in the flat model {£2k = 0). The 
results are shown by the contour lines in the Qm — w plane for D2A (top left panel) 
and D2B (top right panel). As clearly seen in the figures, the relative errors are 
small, less than 4% for 0.1 < flu < 1 and —0.5 < w < —1.5. The errors become 
larger for larger w. This tenancy suggests for larger w that the dark energy has been 
affecting the expansion rate since long time ago, and hence the large differences 
between D2A,b and Df arise at present. 

Figure [2] is the same as Fig. [IJ but for the time varying equation of state 39 )' 40 ) 

w(a) = wq + w a a (1 — a) . (3-4) 

The results are shown in the wo — w a plane with Qm = 0.28 (= 1 — Qx)- As shown in 
the figure, for large w a , the relative differences become large. This is because the dark 
energy term in the hubble expansion H 2 (a), ,!?x a ~ 3( ' 1+W0 '' exp [(3/2)u> a (l — a) 2 ] , 
becomes large for large w a in the past (a < 1). The relative errors are less than 10% 
for —1.5 < Wq < —0.5 and w a < 3. 

§4. Third-order solution 

Similarly, the third-order solution consists of six terms, as shown by 

<5 3 (fe,a) = D 3AA (a)C AA (k) + D' 3AA (a)C' AA (k) + D 3AB (a)C AB (k) 

+D' 3AB (a)C' AB (k) + D 3BA (a)C BA {k) + D 3BB (a)C BB (k), (4-1) 



with 



C AA {k) = ^J fd 3 q a(q, k - q)6 1 (q)A(k - q), (4-2) 

c A A(k) = ^j d3( i k - <i)ti( k - q)Mq), (4-3) 

C AB (k) = ^ J d 3 q a (q,k-q)5 1 (q)B(k-q), (4-4) 

C' AB (k) = lj d 3 Q k - q)S 1 (k - q)B(q), (4-5) 

C BA (k) = ^f d 3 qP(q,k-q)5 1 (q)A(k-q), (4-6) 

C BB (k) = ljd 3 q (3(q, k - q)5 1 (q)B(k - q). (4-7) 
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There are two additional conditions of 

Jg D 3AA + g£>3AB = 2^?' 

\d'zaa + I^sab + ^D 3BA + ^D 3BB = ±Df, (4-8) 

and hence only four terms in Eq. (|4-1|) are independent of each other. The growth 
factors -D3** are determined by the ordinary differential equations with the boundary 
conditions of -D3** — > a 3 in a — > (see Appendix A) . The middle and bottom panels 
in Figs. [T]and [5] are the same as the top panels, but for the relative differences between 
D3** and D\. The results are shown for D^aa (middle left), D$ab (middle right), 
D3BA (bottom left), and D^bb (bottom right). The relative differences are less than 
7% for 0.1 < Qm < 1 and -1.5 < w < -0.5 and less than 20% for -0.5 < w < 0.5 
and w a < 3. 

Our results of the second- and third-order solutions are consistent with the 
previous results of Bernardeau (1994) for the cosmological constant model (w = — 1). 
Although we presented the results for only the density perturbations, one can easily 
obtain the velocity field perturbations by inserting Eqs. (|3T|) and (|4-ip to Eqs. (|2-1|) 
and O^j) . 

§5. One-loop power spectrum 

The one-loop power spectrum is the linear power spectrum with the leading 
correction arising from the second- and third-order density perturbations, 

P(k,a) = <|<5i(fc,a) + 5 2 (k,a) + 5 3 (k,a)\ 2 ) 

= £>?(a)Pn(fc) + P 22 (fc,a) + Pis(k,a), (5-1) 

where P n = (|5i| 2 ), P22 = (N 2 ) and P 13 = (2Re(5i#$))- The first term is the linear 
power spectrum, and the second and third terms are the one-loop corrections. The 
explicit formulae for P22 and P13 are given in Appendix B. 

One usually approximately apply the one-loop power spectrum in the EdS model 
to an arbitrary cosmological model by replacing the scale factor by the linear growth 
factor, 

P Ed s(k,a) = £)?(a)Pn(fc) + Df(a) [P 22 (k) + Pis(k)] ms , (5-2) 

where the second and third terms are the corrections for the EdS model 27 )' 28 ' 1 (see 
also Appendix B). We compare the two power spectra in Eqs. (|5-ip and (|5-2p in 
order to quantitatively demonstrate the validity of the above approximation. We 
use CAMB (Code for Anisotropics in the Microwave Background) 41 -* to calculate the 
linear power spectrum with the cosmological parameters h = 0.701, Qb = 0.0462, 
Q M = 0.279, n s = 0.96 and a 8 = 0.82, consistent with the WMAP 5yr result. 42 ) 

Figure [3] shows the relative differences in ^22(^)5 Pi3(k), P22(k) + Pi3(k) and 
P{k) between Eqs. (|5-ip and (|5-2p at z = 0. The equation of states are (wo,w a ) = 
(-1.2,0), (-1,0), (-1,2), and (-0.8,0). From the top panels, the error is < 5% 
for P13 while <C 1% for P22. On a small scale, these differences are small. In the 
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bottom left panel, the error diverges at k ~ 0.75/i/Mpc because the denominator of 
P22 + -Pi 3 vanishes there. The approximate formula of PEdS predicts a slightly lower 
value than the correct result, because [p22]Eds(> 0) is almost the same as P22 while 
[Pi3]Eds(< 0) is more negative than P13 as shown in the top panels. However, as 
expected, the difference is very small at less than ~ 1% for k < 0.4h/Mpc. Figured] 
is the same as Fig. [3J but at various redshifts of z = 0, 1, 3. Hence, from this figure, 
the EdS model approximation in Eq. (|5 2[) is sufficiently more accurate for higher 
redshifts z > 1. 

Finally, we calculate the shift in the position of the first acoustic peak at k ~ 
0.07/i/Mpc. Dividing P(k) by the no-wiggle model of Eisenstein & Hu (1999), we 
find that the position is shifted by only 0.8%(0.02%) for z = 0(z = 1). 

In this chapter, we calculated the one-loop power spectrum, however it is not 
accurate in the strong nonlinear regime (k > 0.1/t/Mpc). In fact, Jeong & Komatsu 
(2006) found that the one-loop power spectrum coincides with the nonlinear power 
spectrum from the numerical simulation within 1% if A 2 (k) = k s P(k) /2ir 2 < 0.4 
is satisfied. This condition is rewritten as k < 0.12(0.26)/i/Mpc at z = 0(z = 1). 
Hence, in order to extend our analysis to a smaller scale, further analysis of the 
cosmological dependence of the higher-order perturbation theory is necessary. 

§6. Conclusion 

We investigate the third-order density perturbation and the one-loop power spec- 
trum in the dark-energy cosmological model. We present analytical solutions and a 
fitting formula with the general time-varying equation of state for the first time. It 
turns out that the cosmological dependence is very weak, for example, less than 1% 
for k < 0.4/i/Mpc for the power spectrum. However, our results may be useful in 
some cases when one needs a very highly accurate theoretical model of the BAO or 
in the study of the nonlinear evolution on a smaller scale (> 0.4/i/Mpc). 
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Appendix A 

Second- and Third- Order Growth Factors 

The second-order growth factors D-iA,B are determined by the ordinary differ- 
ential equations 



d 2 D 2 
din a 2 a 2 



+ 6 + 



d\nH\ d D 2 



din a J din a a 2 



+ 



n d In H 3 , . 

+ 2-t^ -^M(a) 

am a 2 



D2 
a 2 
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\ da 




fdDA 


I 


V da ) 



+ -f?M(a) 



for D 2B , 



a 



for D 2A , 



with the initial conditions at a = 0: 

D 2 a,b 



_ 1 d D 2 A,B _ 
' da a 2 



(A-l) 
(A-2) 

(A-3) 



For the flat model with the constant equation of state, the solutions are well approx- 
imated as 



D 2A ^ D\ 



D 2B ~ D{ 



1 + \\nQ M \ 



5.54 x 3.40 x 10 



-3 



w 



„ ^ , / 1.384 x 10~ 2 8.50 x 10~ 3 
l + |lnfl M | — + 



(A-4) 
(A-5) 



\ w \ y/\w\ 

within a maximum error of 0.03% for both 0.1 < Q\j < 1 an d —1.5 < w < —0.5. 



Similarly, the third-order growth factors are determined by 



d 2 D 3 
din a 2 a 3 



+ 8 + 



dlniA d L>, 



+ 



ir o dln# 3^ , , 

15 + 3 — -n M (a) 

a In a 2 



18 

y 

15 



dDi 3_ , N Z>! 
da 2 a 



din a J din a a 3 

D 2A ,B , 18 dDi d D 2 A,B 



7 da da a 2 



for D 3 AA,3AB 



dD, 



da 
9 dDi 
2~da~ 



d D 2 a D 2A TDidDi 
a- — + 2 — k 



da a 2 

a d D 2B 2 D2B 
da a 2 a 2 



5 a da 



, for D 3B b 



, for D 3B A 



with the initial conditions at a = 0: 



D 3 



d D 3 



a 3 ' da a 3 

For Qk = with the constant w, the solutions are well fitted by 

'8.21 x 10- 3 5.14 xlO- 3 ' 



D 3 aa ^ Di 
D 3A B ^ D\ 



1 + \lnQ M \ 



w 



(A-6) 
(A-7) 

(A-8) 



(A-9) 



(A-10) 



l + |Infl M | 1 - 6+a4Ia H|fl M | - 7 



H 



9.16 x 10~ 3 8.95 x 1(T 3 
- + - 



\w\ 



(A-ll) 



D 3B A =i D\ 



D3BB — D 1 
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1 + |ln^M| L06_0 - 51n| "' 1 (7.68 x HT 3 M - 1.130 x 1(T 2 



1 + |ln^ M | 



2.641 x 1(T 2 1.582 x 10" 2 



+ 



w 



7 

(A42) 
(A-13) 



within a maximum error of 0.05% for both 0.1 < flu < 1 and —1.5 < w < —0.5. 
The other growth factors D' 3AA and D' 3AB can be obtained using Eq. (|4-8|) . 

Appendix B 

Explicit Expressions of P22 and P± 3 



392vr 2 
1 



Here, we present the explicit expressions of the one-loop correction terms P22 
and P13. From the results in §3, we obtain 

P 22 (k,a) = {\5 2 (k,a)\ 2 ) = Dj A (a)P2AA(k)+2D 2A (a)D 2B (a)P2AB(k)+Dl B (a)P 2 BB(k), 

(B-l) 

with 
P2AA(k) 

P2BB(k) 
P2AB(k) 



roc 

25 fc 4 



d h I ^Pn^Pn^k^ + kl-lkk^j Kk2 + kl _ 2kkifI j , 



kfi + k\ — 2k\\i 



2\ 2 



98vr 2 
5 

196vr 2 



k" 



k 4 



dk-i 



dki 



1 



kfi — k 



1 



dfiPuik^Pu Jk 2 + k\- 2kkxn 



(kfj, - h) (k[i + h - 2km 2 ) 
(k 2 + kf - 2khn) 2 



where [x is the cosine between k and k\. 

Similarly for P13, from the results in §4, we obtain 



P 13 (k,a) = (2Re(<yi(fc,o)^(fc,a))) 

= D 3AA (a)P 3AA (k) + D' 3AA (a)P^ AA (k) + D 3AB (a)P 3AB (k) 

+ D' 3AB (a)P^ AB (k) + D 3BA (a)P 3BA (k) + P 3BB (a)P 3BB (k), (B-2) 



with 

P3AA(k) 

P^AA^k) 
P3AB(k) 
P^ABik) 
P3BA(k) 
PzBB{k) 



roo 

k 3 Pn(k) / drPu(kr) (l + r 2 ) , 
Jo 

I roc 

^k 3 P n (k) J drP n (kr 

1 f°° 
27^ k3Pn ^ J drP n(kr) (1 + r 2 ) , 



54vr 2 
1 

24vr 2 



l + 4r 2_ r 4 + 1^2_ 1 N3 ln 

2r ' 



r + 1 



r — 1 



27vr 2 
1 

168^ 2 ' 
4 



k 3 P u (k) / drPu(kr)r 2 , 



k 3 P u (k) / drP n {kr) 



^( 1 _ 4r 2 _ r 4) + ( r2 



r + 1 



r — 1 
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By setting D2 = a 2 and D3 = a 3 in Eqs. (|B-ip and ()B-2|) . the correction terms 
reduce to the result in the EdS model, a 4 [P22(^) + -Pi3(&)]EdS- 
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. 1. Top panels: the contour lines show the relative differences between Dia,b and D\ at 
present (a — 1) in the Qm — w plane. The flat cosmological model and the constant equation of 
state are assumed. The top left (right) panel is the result for D2A (£>2b). Middle and bottom 
panels: same as top panels, but for the relative differences between D3 and D{. 




Fig. 2. Same as Fig. [T] but for the time-varying equation of state, w(a) = wq + w a a (1 — a). The 
results are shown in the wq — W a plane in the flat cosmological model with Qm = 0.28. 
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Fig. 3. Relative differences of P 2 2(k) (top left), Pi 3 (k) (top right), P 2 2{k) + Pi 3 (fc) (bottom left) 
and P(k) (bottom right) between the correct results and the approximate results denoted by 
[■■•]EdS- The cosmological model is consistent with the WMAP 5yr result. 
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Fig. 4. Same as top panels in Fig[3] but at various redshifts of z — 0, 1, 3. 
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